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An evaluation was made of an approximate method for calculating the 
dose rate resulting from fission product gamma radiation after reactor shut- 
down. Safety analyses for advanced nuclear missions in space have used this 
technique for determining radiologically safe operations. This method 
assumes that the dose rate Di(r, t) , received from a shutdown reactor is 
directly proportional to the residual gamma power level at that time. 

To check the validity of this "power ratio" assumption, energy- 
dependent detector responses were calculated for normalized sources in the 
Perkins and King energy group structure for a SNAP 8 power system on a NASA 
space station. Gamma decay rates were then calculated by using an expanded, 
updated list of isotopic decay data, and from these, actual detector 
responses D 2 (r, t) were found for the SNAP 8 system. The two detector 
responses, Di and D 2 , were then compared at several times after shutdown 
and at detector positions around the space station. An error of several 
thousand percent was found for many detectors at times greater than one 
year; it was also found that an increase in material attenuation produced 
an increase in error and that energy-dependent detector response, such as 
dose rates, showed greater discrepancy between Di and D 2 than did the 
energy flux. 

It was concluded that these discrepancies were caused by a decrease in 
the proportion of high-energy gammas released at longer shutdown times. This 
softening of the decay gamma spectrum with time caused an increased attenua- 
tion of decay gamma energy; invalidating the power ratio assumption. Because 
of the complexity of the space station geometry, these conclusions were 
checked and verified by a point source in infinite water calculation. The 
results indicate that energy-dependent calculations must be made to determine 
decay gamma dose rates for actual reactor configurations. A simplified method 
for making these calculations has been devised. 


This paper summarizes the results from the 
determination of the nuclear environment produced 
by the SNAP-8 reactor on the space station. The 
study was concerned with determination of the gamma 
radiation environment produced by the reactor in 
its shutdown mode, and its variation with respect 
to time elapsed after shutdown. 

The general approach to this and similar shut- 
down reactor radiation problems, which has been in 
use prior to this study, involves two steps. The 
first step is to research the available literature 
and obtain data, generally in graphical form, which 
gives the ratio of residual shutdown gamma power to 
operating power, versus the two parameters of reac- 
tor operating time prior to shutdown and the time 
elapsed after shutdown at which one is interested 
in determining the radiation environment. The 
second step is to assume that the gamma dose rate 
is directly proportional to the residual gamma 
power level, determine the constant of propor- 
tionality, and use it to determine the gamma dose 
rate from the residual gamma power level found in 
the first step. 

The accuracy of this approach depends on the 
accuracy of the graphs from which the gamma power 
ratio is found and the reliability of the assump- 
tion that the gamma dose rate is directly propor- 
tional to the residual gamma power level. Both of 
these matters were investigated. The investigation 
has yielded new power ratio curves which incorporate 


more recent isotopic data than had been used pre- 
viously, a determination of the problems encoun- 
tered, and the error incurred in relating gamma 
dose rate directly to residual gamma power, plus 
a simplified method of performing these calcula- 
tions in the future. 

In addition, a radiation transport analysis 
was made of the nuclear power system for the space 
station, using the Point Kernel technique to deter- 
mine the radiation attenuation characteristics of 
the space station. The results of this portion of 
the study show the relative magnitude of the geo- 
metric and material attenuation, by energy group, 
of gamma radiation in and around the nuclear- 
powered space station. These results will be very 
useful in later investigations in that the gamma 
radiation environment at any of the detector posi- 
tions may be quickly found for a U-235 reactor by 
"de-normalizing" them. i.e. , by multiplying the 
normalized detector response in each energy group 
by the ratio of the actual source in that group to 
the source used in the normalized calculation. 

This process is applicable to future analysis of a 
U-235 reactor at any power level in both the operat- 
ing and post-shutdown phases. 

REVISION OF POST-SHUTDOWN POWER 
VERSUS TIME DATA 

The General Dynamics data on post-shutdown 
gamma power as a function of reactor operating 
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time and time elapsed after shutdown (ref. 1) has 
proven to be a very helpful tool in performing 
radiation analyses on SNAP-type nuclear power 
systems. The original information was compiled 
for a NERVA-type nuclear propulsion reactor, but 
if the same reactor fuel is used (U-235) , the data 
may be applied to any nuclear system. There are 
a few disadvantages in using this material; first, 
it is based on fission product data which was com- 
piled in 1958 and 1959, and, second, the graphs 
contain information for reactor operating times 
only up to 10 4 seconds (approximately 3 hours) . 

It was felt that a major improvement could be 
made to the post-shutdown gamma power data if 
these data were updated to include the latest fis- 
sion product decay information and if results for 
longer reactor operating periods were incorporated. 
References 2 and 3 list the data in use up to 1963, 
with reference 3 having the more recent informa- 
tion, including the fission product data for 125 
nuclides . 

Burrell and Watts (ref. 4) supplied more up- 
to-date fission product data, which included a 
library of decay data for 200 nuclides, plus a 
computer program for utilizing the data. Modifi- 
cations were made to the code, and an output 
plotting technique was written and added to it. 
Figure 1 shows the results of the calculation with 
this code in the form of total shutdown gamma 
power versus time after shutdown and operating 
time. Operating times from 1 to 1 x 10 9 seconds 
(approximately 30 years) were included in the cal- 
culations. The time after shutdown began at 
1 x 10 2 seconds, since this was the minimum that 
could be used from the basic input data without 
extrapolation, and it was felt that periods of 
time shorter than this were not of interest at 
present. 

Figures 2 through 8 display similar data, but 
in this case, separately for the seven Perkins and 
King energy groups. By presenting the data in 
this format, with a separate graph for each energy 
group rather than for each reactor operating time, 
the effect of operating time and shutdown time is 
more apparent. In this form, the results are more 
applicable to future calculations, as described 
in the next section. 


EVALUATION OF POST SHUTDOWN 
DOSE APPROXIMATIONS 

An evaluation was made of an assumption used 
in earlier radiation environment studies of reac- 
tors for nuclear missions. This assumption is that 
the dose rate received from a shutdown reactor is 
directly proportional to the residual gamma power 
level. For example, at a given time after the 
reactor has been shut down when the decay gamma 
power level is P(ti) watts, a calculation deter- 
mines the dose rate at a detector to be Di(r, t\) 
rads. Then, at any other time after shutdown when 
the decay gamma power level is P(t 2 ) watts, the 
dose rate at the same detector can be found from 
the simple equation 

Dl(r ’ tz) = Dt p(t 1 > l) p(t2) 


Note that the ratio of dose rate to residual power 
level at t becomes the constant of proportionality 
by which the new residual power level is multiplied. 
The obvious advantage of this technique for calcu- 
lating decay gamma dose rates is its simplicity. 

Only one transport calculation need be performed; 
it can then be scaled to apply to any time after 
shutdown. 

The following approach was used to check the 
validity of this "power ratio" assumption. As with 
the original assumption, only one transport calcu- 
lation was made for each detector point; however, 
the desired information from this calculation was 
the dose rate as a function of energy. In particu- 
lar, the doses, Di , received from a unit source 
(1 MeV/sec) in each of the seven Perkins and King 
energy groups were calculated. Then, using the 
revised gamma energy release data, F^, as described 
in the previous section, the dose at detector point 
r and at time t was calculated by 

7 

D 2 (r, t) = l Ii(t) Di(r) , 
i=l 

where the summation is over the seven Perkins and 
King energy groups. 

The doses calculated, Di and D 2 , were then 
compared for several times and detector points 
around the space station. The discrepancy between 
doses seemed to depend on three different factors: 
time, position, and type of response. First, the 
error between the two methods increased greatly 
with increasing time since shutdown, reaching over 
1000 percent for some detectors at times greater 
than 1 year. Also, at a fixed time, the error 
became larger as the amount of attenuating material 
between source and detector increased. Finally, the 
error changed at the same detector point and at the 
same time for different response units; for example, 
rads (tissue) versus rads (silicon) versus energy 
flux. Examples of each of these three results are 
given below. 

For the detector located 170.7 centimeters 
above the center of the SNAP reactor, the tissue 
dose rate error between Di and D 2 is shown in 
Table 1. The percent error did not change signi- 
ficantly after 4 x 10 B seconds. 


TABLE 1. TISSUE DOSE RATE 


Time After Shutdown (sec) Percent Error 


1 

X 

10 2 

0 

1 

X 

10 3 

26 

1 

X 

10 4 

66 

1 

X 

10 s 

283 

1 

X 

10 6 

155 

1 

X 

10 7 

590 

1 

X 

10 8 

509 

1 

X 

10 9 

771 

1 

X 

10 10 

770 
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For 1 x 10 5 seconds after shutdown, typical 
detectors showed the discrepancy between D x and D 2 
for dose rates in rads (tissue) per second (table 2). 
These detector positions are also ordered in increas- 
ing amount of attenuating material between source 
and detector. 

For the detector located at 170.7 centimeters 
above the center of the SNAP reactor at 1 x 10 5 
seconds after shutdown, the energy flux showed an 
error of 283 percent, while the tissue dose showed 
an error of 214 percent. Similarly, at 1 x 10 9 
seconds after shutdown, the energy flux showed an 
error of 771 percent, while the tissue dose showed 
an error of 567 percent. 


TABLE 2. DISCREPANCY BETWEEN D x AND D 2 


Detector 

Number 

1 

2 

3 


Relative 
Number of 
Mean Free 
Paths 

1.0 

1.47 

7.90 


Percent 

Error 

214 

333 

1,413 


The physical phenomena causing the large dis- 
crepancies can be understood by examining the gamma 
power versus shutdown time curves (figures 1 through 
8) . These curves show that at early times after 
shutdown most of the gamma energy is coming from 
the higher energy gamma groups. As the time after 
shutdown increases , the percentage of gamma energy 
emitted in the higher energy groups decreases. 

For example, approximately one day after shutdown 
(10 5 seconds), the energy emitted in the highest 
energy group, VII, is only 0.15 percent of the 
total. At 10 2 seconds after shutdown, energy 
group VII contributed 30 percent of the total 
gamma energy release. As the time after shut- 
down approaches 10® seconds (approximately 3 
years), only energy group II contributes signifi- 
cantly to the total decay gamma energy release. 

Group II covers the energy range from 0.4 to 0.9 
MeV. In this energy group, the energy absorption 
coefficients for all materials are greater than 
those for the higher energy groups. Hence, this 
softening of the decay gamma spectrum with 
increasing time, coupled with the greater absorp- 
tion coefficients for softer gammas, caused an 
increased attenuation of the total decay gamma 
energy. For detectors that are shielded from the 
decay gamma source, e.g., in a shutdown SNAP-8 
reactor, the power ratio assumption for calculating 
the shutdown gamma dose would not be expected to be 
very accurate. In fact, the accuracy of the doses 
calculated by this technique should become poorer 
as time increases, because of the softening of the 
decay gamma spectra, and, as the shielding increases, 
because of the selective attenuation of the lower 
energy gammas. Also, these phenomena would cause 
the energy-dependent response functions to show 
variations . 


TABLE 3. TIME DEPENDENT DECAY GAMMA DOSE COMPARISON 


Distance 

from 

Source 

(cm) 



Percent Error at 

Time After Shutdown (sec) 



Response 
Functi on* 

10 ?t 

10 3 

10*’ 

10 s 

10 6 

10 7 

1O 0 

- 

10 10 


EF 

0.0 

0.52 

1.11 

0.91 

0.33 

1.13 

0.82 

2.49 

2.49 


TDR 

0.0 

3.57 

6.53 

7.59 

6.91 

12.10 

11.90 

13.90 

13.90 

10 

EF 

0.0 

0.74 

1.61 

1.62 

0.95 

2.16 

1.84 

3.50 

3.50 

TDR 

0.0 

3.76 

6.97 

8.20 

7.44 

12.90 

12.70 

14,70 

14.70 

25 

EF 

0.0 

0.69 

1.59 

2.83 

2.17 

1.25 

1.32 

0.89 

0.89 

TDR 

0.0 

2.60 

4.37 

4.60 

4.99 

10.60 

10.50 

13.00 

13.00 

50 

EF 

0.0 

7.46 

16.30 

23.30 

18.90 

29.80 

29.20 

26.40 

26,40 

TDR 

0.0 

3.26 

8.04 

12.40 

8.93 

11.40 

10.90 

7.80 

7.79 


EF 

0.0 

29.90 

70.90 

105.00 

87.20 

253.40 

248.80 

268.10 

; 268.10 

100 

TDR i 

0.0 

24.20 

56.60 

83.90 

69.30 

186.10 

182.70 

192.00 

| 191.90 

200 

EF ! 

0.0 

78.00 

242.20 

379.70 j 

297.80 

3426.40 

3210.50 

>10“ 

>10“ 

TDR 

0.0 

72.60 

217.80 

343,10 

270.10 

2904.20 

2727.60 

7885.10 

7883.30 

300 

EF 

0.0 

111.10 

1446.70 

652.50 

495.90 

'9127.90 

8589.90 

>10“ 

>10“ 

TDR 

0.0 

107.80 

420.80 

624.10 

475.60 

[8396.00 

7840.80 

>10“ 

>10“ 

400 

EF 

0.0 

134.00 

673.80 

819.10 

608.00 

>10“ 

>10“ 

>10“ 

>10“ 

TDR 

0.0 

128.80 

614.20 

806.10 

604.30 

>10“ 

>10“ 

>10“ 

>10“ 

500 

EF 

0.0 

143.00 

794.70 

900.80 

670.00 

>10“ 

>10“ 

>10“ 

>10“ 

TDF 

0.0 

142.10 

781.70 

895.30 

666.10 

>10“ 

>10“ 

| >10“ 

>10“ 


*Detector response are energy flux (EF) in {Mev)/(an 2 -sec) and tissue dose rate (TOR) in 
[rads (T)/secj. 


Although the explanation of the cause of the 
discrepancies between Di and D 2 seemed plausible, 
the complexity of the source term and the shield 
geometry in the SNAP-8 reactor system required a 
check of the validity of the above conclusions. 

This check was made by calculating the dose from 
a normalized point source (1 MeV/sec) in each of 
the Perkins and King energy groups in an infinite 
water medium. Detectors were located at nine 
different radii from the point source, beginning 
at 5 centimeters (5 g/cm 2 ) and ending at 500 centi- 
meters (500 g/cm 2 ). Five different responses were 
calculated; however, only the results of two of the 
response calculations will be given. These 
responses are the energy flux (MeV/cm 2 -sec) and 
tissue dose rate [rads (T) /sec] , The effect of 
time after shutdown and shielding on the tissue 
dose rate is very similar to the effect on the 


tResults normalized to 10 2 seconds shutdown time. 
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other three responses. Table 3 shows that, in 
general, the results of the. decay gamma point 
source support the conclusions reached above. 

The percent error increases with increasing time, 
but not monotonically. The error at 10 s seconds 
after shutdown decreases because gamma group IV 
dominates gamma group II at this time. Thus, this 
higher percent of high energy gammas emitted at 
1(T seconds more closely approximates the original 
spectrum (at 10 2 seconds) , hence the decrease in 
percent error between Dx and D 2 . Also, the percent 
error increases for more shielding, i.e., distance 
from the source, except at 25 centimeters. The 
reason for this apparent anomaly is that the build- 
up of dose is greater than the material attenuation 
until some point between 10 and 25 centimeters, 
causing D 2 to be greater than Dx in this range. 

From 25 centimeters on, the original conclusion is 
valid, as expected. 

Considering the original results for the 
SNAP-8 system and the infinite water-point source 
calculations, the following conclusions w erp 
reached: 

• The discrepancies between Dx and D 2 were 
caused by a decrease in the proportion 
of high-energy gammas released at longer 
shutdown times. 

• The softening of the decay gamma spec- 
trum with time causes an increased 
attenuation of decay gamma energy, 
invalidating the power ratio assumption. 

• Energy— dependent calculations need to be 
made to determine the decay gamma dose 
rate at various times after shutdown. 

Obviously, the energy-dependent calculations 
of the dose rate are more time-consuming than using 
the power ratio assumption. However, by using 
detector positions for which a normalized source- 
point kernel calculation has been made, and by 
using figures 2 through 8, the dose rates can be 
calculated as described below. 


Suppose the reactor has operated at a power 
level of P watts for x seconds. The conversion 
factor, C, for data taken from figures 2 through 8 
is computed by 

C " 6.25 x 10 12 [ (MeV/sec) /W] • P (W) * x (sec) 

38 6.25 x 10 12 • P • x MeV. 

For time of shutdown, t, and for each energy group, 
determine the corresponding ordinate values from 
figures 4 through 8. If x, the operating time, 
does not equal the operating times on the figures, 
then use logarithmic interpolation between the 
given values of x. The ordinate values are the 
normalized shutdown gamma power, GP i# for each of 
the seven energy groups. Then, for the selected 
detector position, r, and response, the dose rate 
at shutdown time, t, operating time, x, and power 
level, P, is computed by 

7 

D (r,t) = C • l GP^t) • D^r), 
i=l 


where D^(r) is the energy-dependent detector 
response from the normalized source-point kernel 
calculations. It should be noted that, the above 
process is valid for any shutdown reactor as long 
as the dose rates are determined for a normalized 
source (1 MeV/sec) in each of the Perkins and King 
energy groups. 
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SHUTDOWN GAMMA POWER 










FIGURE 4. - Post shutdown gamma power level energy group (0.9 to 1.35 MeV). 
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TIME AFTER SHUTDOWN (SEC) 

FIGURE 7. - Post shutdown gamma power level energy group (2.2 to 2.6 MeV). 



TIME AFTER SHUTDOWN (SEC) 

FIGURE 8. - Post shutdown gamma power level energy group (2.6 MeV to 00 ). 
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